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Recent developments in the theory of amorphous plasticity point to the central role played by the 
concept of an effective disorder temperature T e //- An athermal dynamics for T e // are proposed in 
the framework of a deformation theory and discussed in light of the recent steady state simulations 
by Haxton and Liu [Phys. Rev. Lett. 99, 195701 (2007)]. The structure of the resulting theory, its 
parameters and transient dynamics are discussed and compared to available data. 



I. INTRODUCTION 

Much recent work has been devoted to the detailed ex- 
perimental, simulational and theoretical analyses of the 
dynamics of low temperature plasticit y in amorphous sys- 
tems p, 1 1 1 s i, 0, s s m mm m m m m- 

The systems of interest include noncrystalline solids well 
below their glass temperature T g , dense granular materi- 
als, and various kinds of soft materials such as foams, 
colloids, and the like. In spite of these efforts there 
remain fundamental open questions that call for fur- 
ther investigation. Some of the most important emerg- 
ing new results pointed to the central role played by 
the concept of an effective disorder temperature T e f / 
0, i, 01, EH, 03 01 d, EE E3- It has been proposed 
that although conventional thermal fluctuations are of lit- 
tle importance for thermal temperatures well below the 
glass transition temperature T g or absent in purely ather- 
mal systems, the state of configurational disorder of the 
deforming system can be characterized by an effective 
disorder temperature T e ff that controls configurational 
fluctuations However, a well-established equation of 
motion for T" e // is still missing. The aim of the present 
work is to discuss such an equation of motion and to 
rationalize its structure based on general considerations. 
The recent steady state simulations by Haxton and Liu 
(HL) are shown to be consistent with the equation 
and are used to determine important ingredients of the 
theory. Predictions for transient dynamics are presented 
and shown to agree with experimental findings. Finally, 
open questions are discussed. 

There has been some recent discussion in the litera- 
ture regarding the mathematical rigor, validity and gen- 
erality of the concept of an effective temperature. It 
turns out that some measures of an effective temperature 
agreed numerically with one another, while others did 
not [Til ]. Furthermore, an exact solution of a simple non- 
equilibrium model shows that there exists a full hierar- 
chy of effective temperatures, that although might agree 
rather well numerically, have different analytic forms 
[l^ ]. Therefore, it seems that for non-equilibrium systems 
there are many ways to define a temperature that quanti- 
fies deviations from thermal equilibrium, many of which 
yield numerically consistent values (for reasons that are 
not yet well understood), but some do not. It is thus nat- 
ural to ask how can one proceed to develop a theory based 



on a concept that is not yet completely well-founded? 

My answer is that there exists ample evidence for the 
existence and importance of an intensive macroscopic 
state variable that characterizes the state of disorder 
of driven amorphous systems (see below) and therefore 
one should try to make progress in a phenomenologi- 
cal way based on general considerations and experimen- 
tal/simulational findings. In saying so, I do not mean 
that we should give up trying to develop a rigorous ef- 
fective temperature formalism. On the contrary, this is a 
major and fundamental challenge. My main point is that 
one can make significant progress even in the absence of 
such an exact formalism. 

From a fundamental point of view it is clear that a 
deformation theory of amorphous systems must include 
a measure of disorder as a basic dynamical ingredient; 
structural disorder is an essential feature of these sys- 
tems. Furthermore, we have a lot of evidence for struc- 
tural evolution during deformation [2(|. We also know 
that deformation history depends on the system prepa- 
ration procedure, which is naturally accounted for by 
assigning different states of initial disorder to different 
preparation procedures [f| 0, 0, H, H(j. In addition, 
there are experiments in which significant structural dif- 
ferences between a material within a shear band and a 
material outside the band were observed, indicating a 
structural sensitivity to the rate of deformation (or the 
rate of energy dissipation) (2(j. In fact, there is evidence 
for structural softening, in addition to thermal softening, 
within shear bands [20, Hl|. Recent experiments show 
that macroscopic quantities like the shear modulus are 
functions of the deformation [22], a point that will be 
elaborated on later. 

An excellent example of these effects was given re- 
cently in a series of simulations by Demkowicz and Ar- 
gon 0, @, 0, Hi ■ These authors were able to identify and 
quantify local structural features of their deforming sim- 
ulated amorphous silicon. They followed the evolution of 
these structural measures and demonstrated beautifully 
the effect of the system preparation procedure on tran- 
sient dynamics as well as the approach to a unique steady 
state (independent of initial conditions) under persistent 
deformation. Similar evidence for the existence of steady 
state of the disorder was found earlier In my opin- 
ion, all this accumulating evidence is pointing towards 
the need to develop a description of the evolution of dis- 
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order in driven amorphous systems and its relation with 
plastic deformation along physically sensible guidelines, 
even prior to the availability of an exact non-equilibrium 
formalism. 

The term "athermal" is used in this work to refer to 
situations in which spontaneous thermal fluctuations are 
incapable of inducing structural rearrangements. There- 
fore, in these situations plastic deformation occurs only 
in response to external driving forces. For example, such 
conditions are relevant for glasses well below their glass 
transition temperature or for granular media and foams 
where thermal fluctuations are practically nonexistent. 
By using this term we by no means imply that thermal 
vibrations or dissipativc mechanisms like friction in gran- 
ular media or viscosity in foams are not important; on the 
contrary these processes are crucial for our discussion be- 
low, where they provide a means to remove energy irre- 
versibly from the deforming system once structural rear- 
rangements induced by the external driving forces take 
place. The ideas developed hereafter are presented in a 
way that makes explicit reference to vibrational motion. 
However, the resulting framework applies equally well to 
systems where thermal vibrations are completely absent. 

In the next section an athermal dynamics for T" e // are 
proposed and discussed in relation to a deformation the- 
ory. The resulting dynamics are compared to the steady 
state Haxton and Liu (HL) data [To| in Sect. IIIII In 
Sect. |lV]we go beyond the steady state analysis to pre- 
dict some transient effects. Section [V] offers a summary 
and some discussion. 



II. ATHERMAL T eff DYNAMICS AND A 
DEFORMATION THEORY 

The basic starting point for the discussion of the effec- 
tive temperature dynamics is the separation of the total 
amount of degrees of freedom of the system under consid- 
eration to fast vibrational and slow configurational ones. 
This separation, known as the "inherent states" formal- 
ism [231 ] . was shown to provide a proper framework to 
describe glassy dynamics; the phase space dynamics are 
such that the system remains in the basin of a single 
minimum in the potential energy landscape (an "inher- 
ent state" or a "configuration" ) for a long time compared 
to the particles vibration timescale until it makes a tran- 
sition to another basin [23|, [24| . Our goal here is to ex- 
tend these ideas to dissipative, driven (out of mechanical 
equilibrium) amorphous systems. 

In the athermal limit considered in this work the tran- 
sitions between the "inherent states" occur only as a re- 
sult of the applied forces and not due to spontaneous 
thermal fluctuations that are assumed to be inefficient at 
low enough temperatures. The response of such systems 
to the application of external driving forces typically con- 
tains a reversible (elastic) and an irreversible (plastic) 
components. The plastic component, usually termed the 
plastic part of the rate of deformation tensor, is denoted 



by D pl . The plastic work density per unit time is SijD?-, 
where Sij is the deviatoric part of the stress tensor. At 
low temperatures, well below the glass transition one, 
non of this energy can be stored reversibly, implying that 
SijDfj is all dissipated fbH ]. 

In light of the separation of degrees of freedom we can 
write the total heat flux Q as 

Q = Qv + Qc, (2.1) 

where v and c stand for vibrational and configurational 
respectively. Generally speaking, Q v is the energy flux 
that is being removed irreversibly from the deforming sys- 
tem to its surroundings. In atomic or molecular systems 
at a finite T, it is being transferred to the heat bath by 
thermal vibrations. In systems where thermal vibrations 
are absent, some other physical mechanism like internal 
friction in granular media or viscosity in foams, is respon- 
sible for irreversibly removing energy from the deforming 
system. Q c is the rate of change of the energy density 
stored in the configurational degrees of freedom. The en- 
ergy dissipation rate equals the total heat flux, therefore 

SijDl] = Q . (2.2) 

Following Langer [l3| we propose that the configurational 
degrees of freedom are characterized by an effective dis- 
order temperature T e ff, possibly different from the bath 
temperature T, such that 

C eff T eff = Qc , (2.3) 

where c e / / is the configurational specific heat of units ks 
per unit volume. Using Eqs. (|2.f p - (|2.4|) we obtain 

CeffT eff = Dg - Q v . (2.4) 

The quantity Q Vl as explained above, represents the rate 
at which energy is removed irreversibly from the deform- 
ing system to the heat bath. In general, a first principles 
calculation that accounts for the detailed processes that 
involve heat exchanges with the bath is well beyond the 
scope of the present work and is anyway an extremely 
difficult task. However, the athermal limit implies some 
simplifications. In that limit, Q v vanishes in the absence 
of plastic deformations, i.e. when Dfj = 0. Moreover, 

Q v is a scalar that is related to dissipation and has the 
dimension of energy density per unit time. These prop- 
erties are naturally accounted for by assuming that Q v 
is proportional to the plastic power density sy.D?-, i.e. 
that 

Qv = s ij -£)fJp(s i j,Te//,r) , (2.5) 

where g(sij,T e ff,T) is a dimensionlcss function. This 
function has a clear physical meaning as it represents the 
fraction of the plastic power density Sij-DfJ that is be- 
ing transformed into regular heat. For example, when it 
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equals unity all the plastic power density SijD^ is con- 
verted into regular heat and non of it is being stored in 
the configurational degrees of freedom. Substituting Eq. 
(f2~5]) into Eq. (|2l| we obtain 

c effTeff = sy-DfJ [1 - g{ Sljl T eff , T)] . (2.6) 

The structure of this equation is consistent with the orig- 
inal equation proposed by Langer [l3[. Following these 
ideas, we can write the total entropy density S as 



and identify 



S — S v -\- S c 



Qc — T e ffS c . 



The last relation implies that 
1 



T, 



eff 



dS c 
dE c 



(2.7) 



(2.8) 



(2.9) 



where E c is the configurational energy density ("inher- 
ent states" energy density). Eq. (|2.9p can serve as a 
definition of the effective temperature (also called the 
configurational temperature) and in fact was used in re- 
lation to Edwards' hypothesis [25| in Ref. [26[, where 
it was shown that a generalized statistical mechanics for 
the "inherent states" can be developed. More important 
to our discussion here are the results of 0, that show 
that Eq. (|2.9[) yields values of T e tf that are consistent 
with standard thermodynamic relations applied to out of 
equilibrium systems. Note also that the above relations 
imply that for out of equilibrium situations, i.e. when 
T e f f^T, we have 



Q = TS V + T e ff S c + TS , 



(2.10) 



which is an expected violation of a standard thermody- 
namic relation. 

The development up to now focused on a quasi- 
thermodynamic interpretation of the effective tempera- 
ture T e ff. In this spirit, it is natural to ask whether gra- 
dients of T e f f drive a configurational heat flux from high 
to low T e ff regions, in analogy with a regular tempera- 
ture behavior. Alternatively, we can ask whether a diffu- 
sion term of the form D e ffV 2 T e f f should be added to the 
right-hand-side of Eq. (|2.6p , where D e f f is an effective 
diffusion coefficient. From the discussion above it is clear 
that if T e ff diffuses, it must be associated with plastic 
flow, i.e. with D pl ^ 0. Therefore, dimensional analysis 
suggests that D e f f cx D pl . The most natural phenomenon 
in which this behavior should be observed is shear band- 
ing. When shear banding occurs, we expect a high T e ff 
region (inside the band, where D pl is large) to coexist 
with a low T e ff region (outside the band, where D pl is 
small), see for example jl7| . Indeed, the shear banding 
simulations of @ suggest that the width of the shear 
band broadens diffusivel y, i .e. as a function of strain (or 
time) to the power 1/2 [23]. This also implies that the 



rate of broadening is proportional to D pl in the band. A 
systematic analysis of such effects is needed in order to 
provide further support for the diffusivity of T e f / . In this 
work we focus on homogeneous situations, such that we 
can omit a possible diffusion term in Eq. (|2.6[) . 

Eq. (|2.6|) becomes useful only if DfJ and g(sij , T e f / , T) 

are known. We first discuss D v \ in the framework 
of the recently proposed Shear- Transformation-Zones 
(STZ) theory [SEllll- Building on early "flow-defect" 
ideas |2j|, this theory views the plastic deformation of 
an amorphous solid as consisting of finite strain, local- 
ized, irreversible rearrangements that take place in re- 
sponse to applied forces. A major theoretical challenge 
is the calculation of the probability of occurrence of these 
these irreversible "events", the so-called shear transfor- 
mations, or more precisely the rate of events which gives 
rise to a macroscopic plastic rate of deformation Dfj. 
The STZ theory culminates in the following expression 
for D% 



D\\ = exp 



Gstz \ fijjsij) 

kllTeff J To 



(2.11) 



which is valid in the absence of extremely fast transients 
[l5| . The validity of this form was demonstrated recently 
in a shear banding simulation Q. 

Let us discuss the physical meaning of this expression. 
The Boltzmann-like factor constitutes the configurational 
part of D p j , a part that is proportional to the probability 
to find a local configurational fluctuation that might en- 
able a plastic event to take place. More explicitly, in light 
of Eq. (|2.9p and in the spirit of Ref. [26[ , we postulate 
that the probability to observe a free energy fluctuation 
SG > Gstz is given by 



P (SG) oc 



f 



■ exp 



SG 



for SG > Gstz 



kBT e ff \ ksT e ff 

(2.12) 

The free energy Gstz represents the existence of a min- 
imal free energy fluctuation below which a shear trans- 
formation cannot occur at all, making the cumulative 
probability to find a potential site for a plastic event pro- 
portional to 



P (6G)d(6G) = exp 



G 



STZ 



knT, 



B^eff 



(2.13) 



as appears in Eq. (|2.1ip . 

The additional multiplicative term fij(sij)/T repre- 
sents that fact that the existence of a potential site for 
a plastic event does not immediately imply the occur- 
rence of an event. There are also dynamic factors that 
determine the actual rate of irreversible events. First, 
the time scale tq determines the basic attempt period 
and was proposed to be of the order of the vibration 
time scale for atomic systems [l^, [H, E3 • For other sys- 
tems, this time scale is the characteristic relaxation time; 
such a time scale always exists, even when thermal fluc- 
tuations are absent. The other factor fij(sij) quantifies 
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the explicit dependence of D pl on the deviatoric stress 
s. The existence of such a term is anticipated on the 
basis of symmetry as D pl is a tensorial quantity. In fact, 
STZ theory offers a more detailed description of the ten- 
sorial function fij(sij), including an additional internal 
state tensor field that accounts for the transition from 
jamming to flow at the yield stress and for orientational 
memory (l2l . ITil . [3fjj . In summary, as described by Eq. 
(|2.11[) , STZ theory provides an expression for the rate of 
irreversible (plastic) events, incorporating both configu- 
rational and dynamic factors. 

The plastic rate of deformation tensor D pl is incorpo- 
rated into an elaso-plastic deformation theory by assum- 
ing that the total rate of deformation tensor D tot can 
be written as a simple sum of linear elastic and plastic 
contributions 



T^tot 

ij 



where efj is the linear elastic strain tensor |31| 



(2.14) 



This 

equation, upon substituting D pl - of Eq. (|2.11[) . together 
with Eq. (|2.6p . constitute a compact deformation the- 
ory. In order to complete the derivation of the theory one 
must specify the functions g(sij, T e ff, T) and fij(sij). As 
a first principles derivation of these functions is not yet 
available, we have to resort to molecular dynamics simu- 
lations. 



roughly the binary collision time. The effective temper- 
ature was evaluated using a relation between the static 
linear response and the variance of the pressure. The 
crucial point is that this measurement of T e f f was shown 
in [ij to yield values that are numerically consistent with 
our definition of Eq. (|2.9p , making HL data most relevant 
for the developed formalism. 

Since our main interest here is in the athermal limit, 
we focus on the low temperatures data and take the 
T ps 0.2To results, where To is estimated from a Vogel- 
Fulcher fit of the equilibrium viscosity, as a representa- 
tive example. All the results arc given in units of To, the 
particle mass, the small particles diameter, the interac- 
tion energy scale and ks = 1- It is important to note that 
as the total rate of deformation D tot was controlled in 
these numerical experiments and only steady states were 
considered, the deviatoric stress s is always larger than 
the yield stress s y , which is defined as 



s„ — lim s 



(3.2) 



where q = ToD tot . 

The steady state deviatoric stress s°° > s y , as well as 
the steady state effective temperature T^*, were mea- 
sured as a function of the normalized imposed total rate 
of deformation ToD tot . A relation between the two was 
found to be of the form 



III. THE HAXTON AND LIU DATA 

An information of the type we are looking for became 
available recently through the remarkable steady state 
simulations by Haxton and Liu (HL) [l(j ■ These authors 
performed molecular dynamics simulations of a purely 
repulsive binary glass, made of equal proportion of small 
and large particles and sheared steadily under a fixed 
area fraction. We write the tensors in the problem as 



f)Pi 





D pl 



D pl 




(3.1) 



and similarly for D tot ; we stress the fact that there is only 
one independent component in each one of them, mak- 
ing the problem effectively scalar. Note also that under 
fixed area the pressure is expected to vary with deforma- 
tion. Modeling direct pressure effects, for example the 
dependence of the yield stress on the pressure, requires 
additional approximations and assumptions that are not 
essential for the present discussion. We therefore keep 
adopting the assumption made above that plastic defor- 
mation is mainly coupled to the deviatoric part of the 
stress and neglect additional pressure effects. 

The ordinary thermal temperature T in the HL simula- 
tions was varied systematically from well above the glass 
temperature to well below and three orders of magnitude 
of imposed total rate of deformation D tot were studied, 
ranging from To-D*°* = 10~ 5 to ToD tot — 10 -2 , where tq is 



so exp 



7^00 



(3.3) 



where G r « 0.0028 and so ps 0.011 arc independent of 
both T (for T < T ) and T D tot . T^ f itself is plotted 
in Fig. [I] as a function of ToD tot . It is observed that 
T^ff approaches a limiting value in the small ToD tot limit 
and exhibits a marked increase for ToD tot > 10 -3 , an 
effect that was not taken into account in previous works 

These results are completely consistent with the T e ff 
dynamics proposed in Eq. (|2.6|) . The steady state rela- 
tion of Eq. (|3.3[) can be recovered, in the framework of 
the proposed theory, by choosing (putting back fc^) 



ff( s y> T e//,T) « ^exp (- - °^ ) 



with 



v • 
(3.4) 

(3.5) 



Substituting Eq. (pT4"|) into Eq. ((23)) we obtain 

Q v = s^Df.^r exp ( -— — ^ — ) for s > s y . (3.6) 
s V k B T ef fJ 

Let us try to interpret these last few results. The struc- 
ture of Q v , i.e. the rate at which energy is being removed 
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FIG. 1: Tg f vs. log 10 (r ( ,D* ot ) in the simulations of Ref. 
[io| |. Data courtesy of T. Haxton and A. Liu. 



irreversibly from the deforming system to its surround- 
ings, is suggestive. The plastic power SijDf- part was 
discussed above. The remaining multiplicative factor, 
which is the function g of Eq. I]3.4[) . represents the frac- 
tion of that power that flows to the heat bath as regular 
heat. It depends both on T e ff and on the deviatoric 
stress through s. The appearance of T e f f in an exponen- 
tial Boltzmann-like factor gives us some confidence that 
Eq. (|3.6[) presents a physically meaningful result and not 
merely a fudge factor, in light of role played by T e f f in the 
proposed theory. This exponential Boltzmann-like factor 
must be related to the distribution of inherent states that 
determines the probability of increasing/decreasing the 
energy stored in the configurational degrees of freedom 
during irreversible plastic rearrangements, though addi- 
tional theoretical work is needed in order to derive it from 
first principles. An important feature of Eq. (|3.4|) is that 
the fraction of s^DfJ that is converted to regular heat 
is an increasing function of T e ff. However, contrary to 
previous works [H, |T3, [l|| [l7| , g also depends explicitly 
on the stress through s. This dependence is such that g 
decreases with increasing s, if we think of T e f f as fixed 
(which is not the case). The fact that g in Eq. (|3.4j) 
depends both on T e f / and s opens the interesting possi- 
bility, as will be shown below, that during transient dy- 
namics energy will be removed temporarily from the con- 
figurational degrees of freedom, resulting in a decrease in 
T e ff. That possibility corresponds to g(sij,T e ff) in Eq. 
(|2.6[) being larger than unity. This effect was demon- 
strated experimentally in [321 ] and will be discussed in 
detail in Sect. [TVl 

In passing we note that Eq. (|2.6[) must be invariant un- 



pi 



Sij and /), , 



dcr the symmetry operations Sij - 
Therefore, the function g(s,T e ff) in Eq. (|3.4|) (which is 
strictly valid for s>s y ) should be smoothly interpolated 



towards s = 0, with a continuous derivative there, for ex- 
ample using So/\/s 2 + s\ with s\ <C s y . In spite of the 
fact that the steady state data for s> s y , corresponding 
to a finite D tot , cannot determine precisely the behavior 
for s < s y we stress the fact that g is an analytic function 
that shows no special behavior near s—>0. In this partic- 
ular relation it is important to note that in interpreting 
their data HL proposed, following the free-volume the- 
ory of [33l | , an equation for T e ff that is consistent with 
the structure of our Eq. (|2.6[) . Specifically, for the rele- 
vant configuration they proposed (using our notation and 
terminology) that 



Q v ra 2s D pl exp ( - 



G r 



kBT e ff 



(3.7) 



which is indeed identical to Eq. (|3.6p for s > s y . How- 
ever, the tensorial generalization of the last relation in- 
volves replacing D pl with D pl , which is defined in com- 
plete analogy with Eq. (|3.5[) . This generalization ex- 
hibits a non-analytic behavior in the limit D pl — >0. This 
kind of non-analyticity was criticized in 34| as being un- 
likely to arise from any first principles analysis of molec- 
ular mechanisms. Wc therefore prefer not to incorporate 
non-analytic behavior in our theory as long as we are not 
forced to. 

The information about the steady state effective tem- 
perature shown in Fig. [T] should be now used to extract 
the only missing piece in our theory, i.e. the tensorial 
function /y in Eq. (12. li p. It has the following structure 
within the STZ theory [S, H| 



fij {Si 



i)=C(s)(- 



(3.8) 



where C(s) is a dimensionless scalar rate function and 
to is an internal state tensor field that represents the 
orientational properties of the plastic events and accounts 
for memory effects and the transition from jamming to 
flow at the yield stress s y [HI, [HI, H(| . 

Let us discuss the physical origin of the tensorial func- 



tion U 



j). The exponential Boltzmann-like factor 



(that depends on T e ff) appearing in Eq. (|2.1ip is propor- 
tional to the average number of available sites for plastic 
rearrangements. The actual rate of plastic rearrange- 
ments, which determines the macroscopic plastic rate of 
deformation D pl , depends also on the current deviatoric 
stress and the recent history of deformation as quantified 
by (s^ , mij ) . The internal state field to quantifies the 
effect of recent history on the deformation and for our 
purposes here should be thought of as a normalized back 
stress. That means that after a plastic event in a certain 
direction took place, further deformation in that direc- 
tion is more difficult, which is represented by the term 



in Eq. Q3.8P [14j |. Therefore, the effective stress at 



a given point is directed according to Sy /s 



This 



expression contains all the tensorial information about 
the symmetry of the deformation, representing its direc- 
tion (but not its magnitude). In an isotropic material, 
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the deformation should follow the symmetry of the lo- 
cal stress Sij] however, recent deformation carries some 
oricntational memory that is represented by rriij whose 
principle axes do not necessarily coincide with those of 

Sij. 

The remaining multiplicative factor in Eq. (|3.8j) . C(s), 
represents the effect of the magnitude of the stress on the 
rate of actual plastic events. This function is certainly 
the most phenomenological part of the STZ theory. A 
first principles derivation of such a material dependent 
function is probably unrealistic and it seems that the best 
approach to study it is by microscopic simulations. In the 
present context we only aim at showing that there exists a 
function C(s), with the general appropriate features, that 
can describe properly the numerical data of HL. First, we 
expect C(s) to be a monotonically increasing function of 
s. Furthermore, following [l6j], we expect that for small 
stresses the dynamics of plastic rearrangements involve 
barrier crossings between different inherent states in the 
energy landscape. Therefore, for s of the order of s y or 
less, we expect activated dynamics to dominate, resulting 
in C(s) that is at least an exponential function of s. For 
higher stresses, we expect the notion of barriers to be 
irrelevant and the rate of plastic rearrangements should 
depend on some other dissipative mechanism, resulting 
in a weaker dependence of C(s) on s. In addition, C(s) 
is a symmetric function that vanishes for s = in the 
athcrmal limit, representing the fact that there cannot be 
irreversible events in a direction opposite to the applied 
force in the absence of efficient thermal fluctuations [3] . 
This property is not relevant for the HL data since they 
report on steady state results with s>s y . 

To prepare the equations for the comparison with the 
HL data, we note that in a simple shear configuration 
the fixed-points of the only non-vanishing independent 
component of m are fl5| 



m (s) 



s/\s\ 

Sy/s 



if H 

if Isl 



< Sy 

> s v 



(3.9) 



In [15| it was shown that m relaxes to one of its fixed- 
points on a time scale much smaller than the typical time 
scale of the T e j f dynamics, enabling us to replace the m- 
dynamics by one of its fixed-point. Therefore, for steady 
state conditions under simple shear, Eq. (|3.8p simplifies 
to 



m (s) 



(3.10) 



We thus conclude the for the purpose of our analysis of 
the HL data, we can use the s > s y branch in Eq. (|3.9[) 
such that Eq. (|2.11[) reduces to 



D pi 



exp 



Gstz \ C(s) 



effj T 



(3.11) 



We now proceed to discuss the estimation of C(s). For 
that aim we define /3 = Gstz /G r and then substitute Eq. 
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FIG. 2: C as a function s°°/s y for /3 = 1 (diamonds), /3 = 1.25 
(squares) and /3= 1.5 (circles). 



(|3~3j) into Eq. (|3~TTj) to obtain 



ToD* 



C = 



1 



(3.12) 



where D tot is considered as a function of s°° . From Eq. 
(|3.12p we learn that C depends both on the yield stress s v 
and p. Estimating s y from the data is a rather subtle is- 
sue as it involves an extrapolation toward the ToD tot — +0 
limit. However, for the sake of our discussion, the exact 
value is not crucial and we choose to estimate s y from Eq. 
(pHj) with T^j(T D tot ->0)~0.001 as the low strain rate 
limit in Fig. [1J The parameter j3, on the other hand, can- 
not be calculated independently, implying that C is not 
uniquely determined by the data. However, we do expect 
it to be 0(1) as it is unlikely that two dramatically differ- 
ent energy scales enter the problem. More interestingly, it 
turns out that a physical constraint bounds f3 from above. 
The idea is that, as explained above, the transitions rate 
function C is a non-decreasing function of the stress s 
since there cannot be less irreversible events in the direc- 
tion of the applied force as its magnitude increases. In 
Fig. [UC is plotted as function of s°°/s y for /3 = 1, 1.25, 1.5 
using Eq. (|3 . 1 2|) . It is observed that for j3 = 1.5 the curve 
becomes non-increasing. Therefore, (3 is constrained to 
be smaller than about 1.5, which seems consistent with 
our original expectation. If indeed f3 ~ 0(1), then the 
rate function C is in the range shown in Fig. [2] Note, 
however, that we cannot derive a strict lower bound for 
/?. Note also that C should be properly interpolated such 
that C— >0 for s— >0, as required in the athermal limit. 

The results presented in Fig. ([2]) follow our general ex- 
pectations regarding the generic features of the material 
dependent function C(s). It is observed that it is a mono- 
tonically increasing function of the stress and that this 
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dependence is stronger for stresses in the vicinity of s y 
than for higher stresses. We conclude that the numerical 
data of HL is completely consistent with the proposed 
theory and tightly constrain those parts of it that at the 
moment cannot be derived directly. In the next section 
we use the resulting theory to predict transient dynam- 
ics that go beyond the steady state data reported in [l(| • 
An important feature of our theoretical structure is that 
the steady state of T e ff is not obtained directly using 
Eq. (|2.6p . which with Eq. (|3.4[) . provides only a relation 
between s°° and T^* . The value of T e °^ is obtained only 
after s°°(T^ f ) is substituted in Eq. ([57TT]) . This implies 
that C(s) in Eq. (|3.11|) is the fundamental quantity that 
determines T£5j. Therefore, additional efforts are needed 
in order to gain a deeper theoretical understanding of this 
function. 
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IV. TRANSIENT DYNAMICS 

The analysis of the previous section was based on 
steady state information alone. In this section we demon- 
strate some generic new transient features of the pro- 
posed T e ff dynamics. For that aim we reduce Eq. (|2.14j) 
to its simple shear form and use the linear elastic relation 



s/2/i, where fi is the shear modulus, to obtain 



s = 2fi [D tot - D pl (s,T t 



eff) 



(4.1) 



with 



TQ 



^( s ,T e// ) = exp - 



Gstz 

ksTeff 



c ( s ) n _ m o( s ) 



(4.2) 

Note that Eq. (|3.10[) was used to allow for arbitrary 
values of s. For simple shear, Eq. (|2.6[) (with Eq. (|3.4p 
and a proper interpolation for \s\ <s y ) reads 



°eff T eff = 



2,D* - 2 : dp% 

V s2 + 4 



exp 



G , 



k B T e ff 



(4.3) 
with s± <^s y . 

For the sake of illustration we keep using the values 
of so, s y and G r corresponding to the HL data (in HL 
units), as well as Gstz = 2G r , /i = 10s y and c e // = l/3. 
Note that the units of c e // are ks per unit volume. In 
HL units, where ks = 1 and the basic volume is an atomic 
one, c e ff is expected to be of order unity. In addition, 
we choose 



C = 0.001|s/s y | 6 /(l 



\S/Sy\ 6 ) 



(4.4) 



which is a simple function with the proper qualitative 
features discussed above. Note that using instead an in- 
terpolated version of one of the curves in Fig. [5] would 
not affect any of the results to follow. We set the ini- 
tial effective temperature to T e ff(t = 0) = 9.5 x 10~ 4 , a 
quantity that, for example, can be controlled by vary- 
ing the rate of cooling by which the glass transition is 



FIG. 3: The deviatoric stress s/s y as a function of the imposed 
strain 7 = / Q D tot (t')dt' . The arrows indicate the direction of 
the applied D tot . 
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FIG. 4: The effective disorder temperature T e // as a function 
of the imposed strain 7 = / D tot {t')dt'. The arrows indicate 
the direction of the applied D tot . 



approached, see |8| for details. We analyzed Eqs. (|4.ip ~ 
(|4.3p for a loading cycle in which ToD tot = 10~ 4 for the 
first half cycle and T$D tot — — 10~ 4 for the second half. 
Figure [3] shows the normalized deviatoric stress s/s y as 

a function of the external strain 7 = J* D tot (t')dt' . The 
figure exhibits a rather generic behavior where there ex- 
ists a stress peak followed by strain softening in the first 
half cycle and a softer response in the second half cycle 
due to memory effects associated with T e f f [l4[ ■ 

More interestingly, Fig. |4]shows the corresponding plot 
for the effective disorder temperature T e f f as a function 
of the imposed strain 7. The qualitative new feature 
observed is that during transients T e ff might exhibit a 
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non-monotonic behavior. During the first half cycle T e ff 
overshoots the steady state value T^j, while during the 
second half cycle it undershoots T£jj. This effect be- 
comes more pronounced as the imposed rate of deforma- 
tion is increased. 

The physical origin of the non-monotonic behavior is 
the dependence of g in Eq. (|3.4p on the stress s, in ad- 
dition to T e ff. If g were a function of T e ff alone, then 
for all cases in which the initial value of T e f f is smaller 
than the steady state one, Eq. (|2.6[) would predict a 
monotonic increase of T e ff, irrespective of the dynam- 
ics of the stress s. However, the dependence of g in Eq. 
p.4[) on the stress s, allows for a non-monotonic behav- 
ior of T e ff during transients. For example, during the 
first half of the loading cycle discussed above, g(s,T e ff) 
becomes unity while the stress is still decreasing in the 
strain softening part observed in Fig. [3l That implies 
that g(s,T e ff) in Eq. (|3 .4[) temporarily overshoots unity 
such that T e f f becomes negative. This behavior suggests 
a somewhat surprising physical phenomenon in which 
during deformation transients energy stored in the con- 
figurational degrees of freedom is released, i.e. Q c < 0, 
and converted into regular heat such that the amount of 
heat being removed from the system to its surroundings 
is larger than the plastic power SyDfJ. 

The predicted non-monotonic behavior of T e ff during 
transients can be tested directly in laboratory experi- 
ments. For that aim, consider the heat equation for the 
regular temperature T 

cf = g(s, T e// )s y Z»fJ + DV 2 T , (4.5) 

where D is the thermal diffusion coefficient and c is the 
thermal specific heat of units ks per unit volume. Note 
that thermoelastic effects were neglected here. Under 
sufficiently fast deformation the heat conduction term 
can be neglected, leading to adiabatic conditions. There- 
fore, a direct time-dependent measurement of both T and 
SijD?j allows the determination of g(s,T e f f). In fact, 
such a procedure was realized and reported in [32| . In 
this work the temperature rise during high strain rates 
deformation of a glassy polymer was measured. The real- 
time temperature measurements were obtained by using 
embedded thermocouples and transient high strain rates 
of up to 8000s" 1 were obtained by impact boundary con- 
ditions. The plastic power was evaluated by subtracting 
the elastic power from the total one. Then, the time- 
dependent temperature record and the known value of 
c in Eq. (14. 5|) were used to determine g (denoted by 
Pdiff m [32j|V The most striking result was that g be- 
came larger than unity during fast transients. Moreover, 
the effect took place during the strain softening part of 
the stress-strain curve and was more pronounced as the 
strain rate was increased, all in perfect agreement with 
the predictions of the proposed theory. Therefore, the 
results reported in [32| clearly demonstrate an important 
physical effect that emerges from our equations. 

The discussion that follows Eq. (|4.5|) points to a pos- 
sible limitation of the HL simulations. The idea is that if 



strain rates of ~ 10 4 s~ 1 already lead to adiabatic condi- 
tions and therefore to a non-negligible temperature rise, 
then the enormous rates of deformation reached in HL 
simulations (up to 1% per unit atomic vibration, see Fig. 
[IJ must lead to melting of their glass due to the inabil- 
ity of the material to conduct heat sufficiently fast. This 
physical effect was most probably avoided in the simula- 
tions by allowing the numerical thermostat to pump en- 
ergy from the deforming system at whatever rate needed 
in order to keep T fixed, effectively using an unrealis- 
tic thermal diffusion coefficient. Another manifestation 
of the same problem can be seen in Fig. [5] where the 
stress can reach values more than an order of magnitude 
larger than the yield s y . This observation may ques- 
tion the relevance of the high rates of deformation part 
of the HL data to understanding plasticity of real ma- 
terials. However, even if we neglect the very high rates 
of deformation (and very high stresses) part of the HL 
data, the remaining part still provides an excellent test 
for the ideas presented in this work and can be used to 
isolate novel and important physical aspects of plasticity 
of amorphous systems. The novel non-monotonic tran- 
sient effects, which are supported by the available exper- 
imental results of [32| is one of them. 

It is important to stress that the non-monotonic be- 
havior of T e ff during transients might be very impor- 
tant for understanding strain localization. In [I7| it 
was shown that shear banding initiation takes place only 
during transients, while the steady states are linearly 
stable. However, the authors of [17| used an equation 
for T e ff that precludes a non-monotonic behavior and 
a strain rate dependent T^j. Redoing the calculation 
with the T e ff dynamics proposed in this work, i.e. Eq. 
(|4.3p . might shed some new light on this important phe- 
nomenon. Moreover, using Eq. (|4.3[) the authors of [36[ 
found a linear instability in which an expanding circular 
cavity becomes unstable against the formation of crack- 
like localized fingers, a result that might be relevant for 
understanding "brittle" versus "ductile" behaviors. This 
clean linear instability was not found for a T e f / equation 
that allows only for a monotonic approach to steady state 
and see also 133 for details. 



V. SUMMARY AND DISCUSSION 

In this work an attempt to formulate and rationalize an 
equation of motion for an effective disorder temperature 
T e f f was made, building on earlier ideas of Langer [l3T | . 
The resulting equation, coupled to the predictions of the 
STZ theory [IJ|, suggests a compact deformation the- 
ory for athermal amorphous systems. One of the results 
obtained is that this theoretical structure is consistent 
with the steady state simulations of for low tem- 
peratures. However, the functions g(sij,T e ff,T) in Eq. 
(|2.6|) and C(s) in Eq. (|3.8|) were not derived from first 
principles, but were extracted from the numerical data in- 
stead. Understanding how to derive these functions from 
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a more fundamental point of view is a major challenge. 
The present work focused on athermal conditions; un- 
derstanding how to treat theoretically the contributions 
of ordinary thermal fluctuations and the interaction of 
the deforming system with the heat bath for higher tem- 
peratures T is yet another challenge. Some steps in this 
direction were made in [l6l . [3?| . 

Several recent works proposed alternative T e f f dynam- 
ics. We have already mentioned that an equation very 
similar to the one proposed here was suggested in [10l ] . 
Moreover, the equations used in [n| [35[ are fully consis- 
tent with the general structure of Eq. (|2.6[) . but with 



g(s iv T ef} ) 



T, 



eff 



T™ f (T DP l ) 



(5.1) 



instead of our Eq. Here T^ f {r D pl ) is the function 

plotted in Fig. [TJ In fact, in [16[ it was shown that this 
choice is consistent with the very same HL data analyzed 
here, implying that neither the choice of Eq. (|3.4p nor 
that of Eq. (|5.1[) is unique. However, our transient dy- 
namics analysis provides a way to distinguish between the 
two approaches. Eq. (|5.1[) immediately implies that T e // 
reaches its steady state value monotonically, whereas we 
have demonstrated explicitly that Eq. (|3.4[) may lead to 
non-monotonic transients that were observed experimen- 
tally in [32I ]. This difference, as explained above, may 
be very important for understanding strain localization 
that was shown to occur during transients [l6|, [3(| • More 
importantly, there is a fundamental difference between 
the two approaches. When using Eq. (|5.ip . the function 
T£j j{tqD p1 ) is an input to the theory and is not calcu- 
lated within the theory. However, in our approach, using 
Eq. (|3.4p , T£f f (T D pl ) is determined by the theory once 

C(s) is known, making the stress dependent rate function 
C(s) a fundamental quantity. Therefore, we still have to 
understand whether the function ^{tqD p1 ) is a funda- 
mental and general quantity, depending only on D pl , but 
not on the details of the rearrangements, as was claimed 
in [l|| ; or that the stress dependent rate function C(s) 
that captures much of the microscopic physics of the irre- 
versible rearrangements is more fundamental, as implied 
by our approach. 

The present work provides further support for the role 
of the concept of an effective disorder temperature in un- 
derstanding the irreversible deformation of amorphous 



systems. In light of these results, based heavily on com- 
puter molecular dynamics simulations, one must ask how 
these ideas are to be tested in laboratory experiments. 
This raises the question of how to measure T e f f in real ex- 
periments. In fact, a beautiful example for such a labora- 
tory measurement was give in [38j where the mobility and 
self-diffusion of tracer particles within a dense granular 
flow were used to evaluate T e ff through an Einstein re- 
lation. This measurement is directly related to the work 
presented here since in (2?J it was shown explicitly that 
such an Einstein relation between the self-diffusion coef- 
ficient and the mobility coefficient yields a T e f f that is 
consistent with our definition of T e ff in Eq. (|2.9p . Other 
possibilities for laboratory measurements of T e f f include 
coupling the system to a low- frequency harmonic oscilla- 
tor and measuring its fluctuations [391 ] or measuring the 
barrier crossing rate of a two-level system weakly coupled 
to the system of interest, as demonstrated recently in a 
simulation [ill ]. 

An entirely different and new possibility for measur- 
ing T e ff is related to the quasi-thermodynamic role of 
T e ff, as stressed in [i| [l5j . In these works T e // was 
used as a state variable in equations of state for different 
physical quantities. For example, in fl5| an equation of 
state relating the mass density p to T e ff was proposed 
for simulated amorphous silicon. Once such equations of 
state become theoretically available one can obtain sim- 
ple "T e f /-thermometers" by measuring standard macro- 
scopic quantities as a function of the deformation and 
relating them to T e f f through the equation of state • 
An indication for this possibility was given recently in 
[22| where the shear modulus fx was shown to be a func- 
tion of the rate of deformation. This quantity can be in 
principle related to T^ f (D pl ) once an equation of state 
of the form fj,(T e ff) is known. This is yet another future 
line of investigation emerging from the present work. 
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